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ABSTRACT 

We present a model for the formation of high-mass close binary systems in the 
context of forming massive stars through gas accretion in the centres of stellar clusters. 
A low-mass wide binary evolves under mass accretion towards a high-mass close binary, 
attaining system masses of order 30-50 M at separations of order 1 AU. The resulting 
high frequency of binary systems with two massive components is in agreement with 
observations. These systems are typically highly eccentric and may evolve to have 
periastron separations less than their stellar radii. Mergers of these binary systems are 
therefore likely and can lead to the formation of the most massive stars, circumventing 
the problem of radiation pressure stopping the accretion. The stellar density required 
to induce binary mergers is 
collisions. 
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1 INTRODUCTION 

The formation of high-mass stars is a large unknown in mod- 
ern astronomy. While our understanding of the formation 
of low-mass stars has improved dramatically over the past 
decade (Larson 2003), we still do not know whether massive 
star formation is basically a scaled-up version of low-mass 
star formation or if it results from a dramatically differ- 
ent process (Stahler, Palla & Ho 2000; Bally & Zinnecker 
2005). The basic questions concern how mass is added and 
how it interacts with the large radiation pressure from the 
high stellar luminosity once the star reaches masses in ex- 
cess of IOMq. Models of scaled-up low-mass star formation 
have repeatedly shown that if the gas contains typical inter- 
stellar dust, then the radiation pressure deposits sufficient 
momentum into the dust (well coupled to the gas) to repel 
the infalling gas (Yorke & Kriigel 1977; Wolfire & Casinelli 
1987; Beech & Mitalas 1994; Edgar & Clarke 2004). 

Potential solutions to this problem include accretion 
through a disc (Yorke & Sonnhalter 2002), accretion over- 
powering the radiation pressure (McKee & Tan 2003), and 
massive star formation through mergers of (dust free) stars 
(Bonnell, Bate & Zinnecker 1998; Bonnell & Bate 2002). 
Disc accretion can help in two ways. Firstly, the matter ac- 



E-mail: iabl@st-and.ac.uk 



cretes through a small solid angle and is thus less exposed to 
the radiation pressure. Secondly, as young stars are generally 
thought to be rapidly rotating, the star can be significantly 
cooler at the equator than at its poles. Thus, there is less 
radiation pressure in the equatorial regions to oppose the 
accretion. Using this last factor in their models as well as 
a sophisticated treatment of the radiation transfer, Yorke 
& Sonnhalter (2002) showed that it is possible to accrete 
up to masses of the order 30 Mq before the radiation pres- 
sure repels the infalling envelope. McKee & Tan (2003) have 
speculated that a very concentrated centrally condensed core 
will provide large accretion rates onto a forming massive star 
which will then overpower the radiation pressure of the grow- 
ing star. There are two potential limitations to this process. 
Firstly, such a core, even if centrally condensed and sup- 
ported by turbulence, is very likely to fragment and form 
a small stellar cluster (Dobbs, Bonnell & Clark 2005). Sec- 
ondly, the implied accretion rates are similar to the models 
of Yorke & Sonnhalter (2002) which show that the infalling 
envelope is indeed repelled once the star attains masses of 
order 30 M . 

The third potential solution involves the merger of 
lower-mass stars in a dense stellar cluster (Bonnell et 
al.1998). While certainly the most exotic of the three, it 
does have the attraction that any dust will be destroyed 
in the lower-mass stars and thus there is no problem with 
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the radiation pressure from the forming massive star. The 
obvious difficulties with this model is that it requires very 
high stellar densities (of order 10 8 stars pc~ 3 ) in order for 
stellar collisions to be sufficiently common to form massive 
stars in less than 10 6 years. Although such densities are 10 3 
times higher than generally found in the cores of young stel- 
lar clusters, they are not inconceivable as the dynamics of 
accretion can induce the core of a cluster to contract signif- 
icantly (Bonnell et al.1998; Bonnell & Bate 2002). 

One inescapable feature of massive stars is that they 
form in rich stellar clusters (Clarke et al. 2000; Lada & Lada 
2003). Even apparently isolated massive stars are best ex- 
plained as being runaways from stellar clusters (de Wit et 
al.2005; Clarke & Pringle 1992). These systems have a reg- 
ular field-star IMF and thus we need to put massive star 
formation into the context of forming many more low-mass 
stars. Models for the formation of stellar clusters, neglect- 
ing any feedback processes, show that the fragmentation of 
a turbulent molecular cloud can form hundreds of stars and 
that they follow a field-star like IMF (Bonnell, Bate & Vine 
2003). Intriguingly, the initial masses of these stars, even 
those that end up being high-mass, are all initially close to 
the Jeans mass of the cloud (Bonnell, Vine & Bate 2004). 
The development of high-mass stars in these systems oc- 
curs due to subsequent accretion onto the forming cluster, 
where the stars near the centre of the potential accrete more 
rapidly and thus attain much higher masses (Bonnell et 
al. 2001a; Bonnell et al.2004). This effect is increased once 
a significant disparity in masses is achieved. It is this com- 
petitive accretion which explains the origin of the IMF while 
maintaining a low median stellar mass in the cluster (Bon- 
nell et al. 2001b;, Bonnell et al.2003). Massive star forma- 
tion in the context of low-mass star formation can thus be 
explained as being primarily due to the environment of the 
forming massive star (Bonnell et al.2004), linking the forma- 
tion of a massive star to the formation of a stellar cluster. 
However we still need to understand how the mass accretion 
overcomes the radiation pressure to form the most massive 
stars. 

A further complication to any model of massive star 
formation is that most massive stars are in binary systems 
(Mason et al.1998; Preibish et al.1999; Garcia & Mermilliod 
2001). These systems are often very close, with separations 
<1 AU, and companions that are also high-mass stars (Gar- 
many, Conti Massey 1980; Bonanos & Stanek 2005). It is 
these close massive binary systems that motivate the present 
study. In Section 2 we explain the difficulties of forming close 
binary systems. In Section 3 we discuss the numerical simu- 
lations of cluster formation. Section 4 explains the formation 
of massive binary systems in stellar clusters and their orbital 
evolution. The potential for binary mergers is discussed in 
Section 5. We discuss the implications of these results in § 6 
and finally, our conclusions are given in Section 7. 



2 THE PROBLEM: FORMING CLOSE 
BINARY STARS 

Forming close binary stars systems is difficult even amongst 
lower-mass stars, If the binary components form through 
fragmentation (eg, Boss 1986; Bonnell 1999), then the Jeans 



radius at the point of fragmentation must be smaller than 
the binary separation, 

i?bin>2 J R J ocT 1/ V 1/2 - (1) 
This implies a high gas density and thus a low Jeans mass, 
M, ~ M.j oc T 3/ V 1/2 - (2) 

This results in the mass of the individual stars being directly 
related to their separation, 

■Rbin OC — , (3) 

such that close systems have very low masses (Boss 1986). 
For example, if the typical 30 AU binary has solar mass com- 
ponents, then a 1/3 AU binary should have components of 
0.01 Mq (Bonnell & Bate 1994). Forming close binary stars 
in situ is therefore difficult as it requires subsequent accre- 
tion to reach stellar masses (Bate 2000). An alternative is 
that the components form at greater separation and then 
are brought together. Recent simulations of low-mass star 
formation in a cluster environment have shown that close 
binaries can result from the induced evolution of wider sys- 
tems (Bate, Bonnell & Bromm 2003a). The binaries evolve 
due to gas accretion, angular momentum lost to circumbi- 
nary discs (e.g. Pringle 1991;Artymowicz et al.1991) and 
dynamical interactions with other stars. Can the same pro- 
cesses explain high-mass close binary systems? 

2.1 Accretion and Binary Evolution 

Accretion onto binary systems has the potential of form- 
ing close systems out of wider systems at the same time as 
forming higher-mass components. In order to see this, let us 
consider the angular momentum of a binary system, 



(4) 



If the accreted material has zero net angular momentum, as 
is expected if it infalls spherically symmetric, then L ss con- 
stant, and the binary separation should be a strong function 
of the mass, 



R bin oc M hi ' E 



(5) 



If instead, the accreted material has constant specific angu- 
lar momentum, the same as the initial binary, then the total 
angular momentum will scale with the mass of the binary, 
L oc M and thus the orbital separation will scale with the 
mass as 



Rbin OC M bi l 



(6) 



In a turbulent medium, as expected here, the angular mo- 
mentum of each parcel of infalling gas is essentially randomly 
oriented. The net angular momentum will then do a random 
walk with increasing mass such that 



£°cM b Y n 2 , 
and thus 

R hin oc M^. 



(7) 



(8) 



Under these basic assumptions, we can see that accretion 
onto a binary system can significantly decrease its separa- 
tion at the same time as it increases its mass (Bate & Bon- 
nell 1997; Bate 2000). 
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3 CALCULATIONS 

The numerical simulations discussed here were first reported 
in Bonnell, Bate & Vine (2003). The simulations used the 
Smoothed Particle Hydrodynamics (SPH) method (Mon- 
aghan 1992). The code has variable smoothing lengths in 
time and in space and solves for the self-gravity of the gas 
and stars using a tree-code (Benz et al. 1991). The initial 
conditions consisted of 1000 Mq of gas in a 0.5 pc radius 
uniform density sphere. The cloud also contains significant 
turbulent motions such that the kinetic energy is equal to 
the magnitude of the potential energy. The gas is isother- 
mal at 10K which implies an initial Jeans mass of IMq. 
Star formation is modeled by the inclusion of sink-particles 
(Bate, Bonnell & Price 1995) that interact only through self- 
gravity and through gas accretion. Sink-particle creation oc- 
curs when dense clumps of gas have p > 1.5 x 10~ 15 g cm -3 , 
are self-gravitating, and are contained in a region such that 
the SPH smoothing lengths are smaller than the 'sink ra- 
dius' of 200 AU. This ensures that the initial sink-particles 
can have masses as low as our resolution limit of Q.IMq. 
Gas particles are accreted if they fall within a sink-radius 
(200 AU) of a sink-particle and are bound to it. In the 
case of overlapping sink-radii, the gas particle is accreted 
by the sink-particle to which it is most bound. The grav- 
itational forces between sink-particles is smoothed within 
160 AU using the SPH kernel. Thus any binary separation 
within 160 AU is an overestimate of the true separation. 
This implies that the frequency of dynamical interactions in 
the simulation is artificially higher than would be the case 
if the binary was allowed to evolve to smaller separations. 
This is offset to some degree by the corresponding decrease 
in the strength of each interaction due to the gravitational 
softening. The simulations were carried out on the United 
Kingdom's Astrophysical Fluids Facility (UKAFF), a 128 
CPU SGI Origin 3000 supercomputer. 



3.1 Reconstructing the binary's orbital 
parameters 

In spite of the fact the gravitational forces are smoothed at 
distances less than 160 AU, we can still hope to extract 
information as to what the true binary separations should 
be. This is possible as the orbital parameters depend solely 
on the system's mass, angular momentum and total energy, 
which are all directly calculable from the simulation. Thus 
we can estimate what the true binary semi- major axis, i? S omi, 
would be in the absence of any gravitational smoothing as 



itational forces. We can then estimate the semi-major axis 
of the orbit as well as the periastron separation, 



Rsa 



GM bil 



(9) 



where J is the specific angular momentum, M is the total 
mass of the binary system and G is the gravitational con- 
stant. We can likewise calculate the eccentricity, e, of the 
orbit from 



1 + 



2EJ 2 

bin 



(10) 



where E is the total specific energy (energy per unit mass) 
of the binary system. This term includes the smoothed grav- 
itational potential that corresponds to the smoothed grav- 



Rt, 



-Rsemi(l c), 



(11) 



of the binary and determine its evolution. Furthermore, we 
can compare the periastron separation to the stellar radii to 
estimate if the binary system would merge to form a more 
massive star. 



4 FORMATION OF HIGH-MASS BINARY 
STARS 

The birth of high-mass stars in the simulation of Bonnell et 
al.(2003) is due to competitive accretion in a cluster envi- 
ronment. The simulation forms numerous sub-clusters which 
eventually merge to form one large system. The high-mass 
stars form in the centres of the dense sub-clusters due to the 
combined potential that funnels the gas down to the centre 
of the system there to be accreted by the growing proto- 
massive star. The clusters themselves form as small-N group- 
ings and grow by accreting stars and infalling gas. Thus, the 
system is initially small and 3-body capture occurs readily 
(Binney & Tremaine 1987). In small-N systems a central bi- 
nary forms when three initially unbound stars pass so close 
that there is violent exchange of energy, with one star being 
ejected at high speed and the other two becoming bound 
to one another. In larger systems this does not work as the 
larger velocity dispersion and smoother gravitational poten- 
tial drastically reduces the probability of having a third star 
sufficiently close to extract the excess kinetic energy. As the 
high-mass stars all form in the centre of the individual sub- 
clusters, they generally have undergone three-body capture 
early in the cluster's growth and are therefore in binary sys- 
tems. 

Once the binary has formed, continuing accretion onto 
it increases the masses of the individual stars. Of equal im- 
portance is the effect of the accretion on the binary's separa- 
tion. As the infalling gas has no correlation with the binary, 
its specific angular momentum is uncorrelated with that of 
the binary's. Thus, accretion does not significantly increase 
the binary's angular momentum and the separation of the 
binary decreases as mass is added (see also Bate 2000) . This 
process is illustrated in Figure 1 which shows the evolution 
of the binary mass versus orbital separation for one of the 
high-mass binary systems formed in the simulation (final 
masses 17 and 10 Mq). The system originates as a low-mass 
wide binary and evolves towards higher masses and smaller 
separations. Once the separations are < 160 AU, the gravi- 
tational smoothing stops the binary from evolving to smaller 
separations. 

The early stages of the evolution are well parameterised 
by i?bin oc M^? as expected for accretion from a turbulent 
medium where each gas parcel has randomly oriented angu- 
lar momentum (see §2.1 above). Unfortunately, the binary's 
orbital evolution quickly enters the regime where the gravi- 
tational forces are smoothed (160 AU). From this point on, 
we use the evolution of the energy and angular momentum 
to determine the evolution of the binary's orbit. 
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Figure 1. The evolution of the binary separation is plotted 
against system mass for one of the high-mass binaries formed in 
the simulation. The system forms as a low-mass wide binary due 
to three-body capture and evolves to a high-mass closer binary 
due to gas accretion. The separation is an overestimate of the ac- 
tual separation due to the gravitational smoothing at separations 
less than 160 AU. 
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Figure 2. The binary's system mass is plotted against the de- 
duced semi-major axis (dashed line) and periastron separation 
(solid line) for the same high-mass binary as in Figure 1. For 
comparison, the separation due to the gravitational smoothing is 
plotted as the dotted line. Note that the evolution is from low- 
mass wide system to a high-mass close system with i? scm i J$ 10 
AU and fl peri < 1 AU. 



4.1 Evolution of the binary parameters 

A binary's orbital parameters depend solely on the angular 
momentum, energy and mass of the binary system (See §3 
above). Thus, as long as we can calculate these quantities, 
we can extract the binary properties and their evolution. 
Figure 2 shows the evolution of the binary's system mass as 
a function of its semi-major axis and periastron separation, 
for the same binary as shown in Figure 1. The evolution of 
the smoothed separation is also shown for comparison. The 
semi-major axis of the binary is generally much smaller than 
the smoothed separation. Thus, from separations of order 
100 AU at masses of several solar masses, the binary evolves 
to a semi-major axis of order 1 AU by the time the system 
has accreted up to 30 Mq. This once again roughly agrees 
with an iibin <x M^? evolution as expected for accretion 
from a turbulent medium. 

The evolution of the periastron separation is also plot- 
ted in Figure 2. from this figure we can see that the binary 
is generally in a highly eccentric orbit where the periastron 
separation is an order of magnitude smaller than the semi- 
major axis, corresponding to an eccentricity of e « 0.9. At 
one point the periastron separation appears to be larger than 
the semi-major axis and this denotes a binary system which 
is unbound. At other points the evolution decreases in sys- 
tem mass indicating an exchange has occurred with a third 
star, or that this third star temporarily passes closer to the 
primary. There is obviously some ambiguity in the evolution 
due to the gravitational smoothing imposed which maintains 
an artificially large cross section for further dynamical in- 
teractions in addition to smoothing out these interactions. 
Nevertheless, we see that the periastron separation decreases 
more dramatically than does the semi-major axis and attains 



separations of order the size of high-mass stars (7? por i < 0.05 
AU). In fact, at several points, the deduced semi-major axis 
is smaller than the stellar radii of the stars. Such systems can 
reasonably be expected to undergo stellar mergers or in the 
very least, tidal forces which will reduce their semi-major 
axes to be of order their periastron separations. 



5 BINARY SYSTEMS AND STELLAR 
MERGERS 

The above evolution of one of the high-mass binaries is a 
typical example of the evolution of the orbital parameters 
under mass accretion. The distribution of binary separation 
as a function of mass is shown for all the high-mass binary 
systems (where at least one star has m > 5Mq) in Fig- 
ure 3. This figure plots the binary's semi-major axis against 
the total binary mass (filled pentagons) and the individual 
component masses against the periastron separations (open 
triangles), linked by solid lines for each system. The more 
massive stars are generally in binaries with semi-major axes 
less than 10 AU and periastron separations less than 1 AU. 
In fact, of the 10 stars at the end of the simulation with 
masses m > IOMq, all ten are in binary systems J? SC mi < 100 
AU and 8 of them have i? SC mi < 10 AU. For the 18 stars 
between 5 and 10 Mq, the numbers are 7 with i? SC mi < 100 
AU and 5 with i? scmi < 10 AU. 

In addition to the large fraction of close binary systems 
amongst the high-mass stars, it is also worth noting that 
generally both components of the binary are high-mass stars. 
Both of these properties are in agreement with observations 
of massive binary systems (Mason et al.1998; Garcia & Mer- 
milliod 2001; Garmony et al.1980). Of even greater potential 



Massive binary stars 5 



3 



r. 




mass Accretion Mass 



Figure 3. The binary semi-major axis is plotted against binary 
mass (filled pentagons) for all the systems containing at least 
one star with m > 5Mq. The periastron separations (open tri- 
angles) are also plotted for these systems against the individual 
component masses, and are joined to their total binary masseand 
semi-major axis by solid lines. Of note is that some binary sys- 
tems have only one individual component plotted indicating that 
the other star has another, closer companion, The dashed line in- 
dicates an approximation of the stellar radii and thus periastron 
separations below this line would force the binary to merge. 



Figure 4. The mass gained in binary mergers, based on hav- 
ing periastron separations smaller than one half of the primary's 
radius, is plotted against the mass acquired through gas accre- 
tion. The final products are plotted as filled pentagons while 
open squares denote merger products which subsequently merged 
with a more massive star. The mergers are estimated by calculat- 
ing the binary properties throughout the simulation and ensuring 
that each star can only merge once with a more massive star. The 
solid line denotes equal mass being acquired through mergers and 
accretion. 



significance is that the periastron separations of the systems 
are dangerously close to the estimated size of the stellar radii 
(assuming a main-sequence mass-radius relationship), 



Rg; 



( M, 
I Mo 



(12) 



Thus, many of these systems could be undergoing mergers, 
while at least one system has 7? pe ri < R*. In fact, this is 
only an instantaneous picture of the binary systems and as 
we can see from Figure 2, the evolution is somewhat chaotic 
and frequently perturbed by passing stars. 

In order to quantify the frequency and relevance of bi- 
nary mergers in high-mass star formation, we have evalu- 
ated throughout the simulation when a binary system has 
a periastron passage closer than half its stellar radii. Us- 
ing this and the condition that each star can only merge 
with a more massive star once (only the most massive pre- 
merger star continues as an independent entity), we have 
estimated the potential for mass growth via binary mergers. 
This merger- acquired mass is plotted in Figure 4 against the 
mass acquired through gas accretion. The stars that sur- 
vive the merger process are plotted as filled squares while 
those that merged during the evolution are plotted as open 
squares. We note from this figure that while most of the 
more massive stars should undergo mergers, these mergers 
are generally with lower-mass stars and thus would not sig- 
nificantly increase the stellar masses. This is not the case 
for at least three of the stars which would more than double 
their masses through mergers. In particular, the most mas- 



sive star would increase its mass from w 30Mq to nearly 
75M Q through multiple binary mergers. Thus, in some cases 
binary mergers could be a significant factor in forming mas- 
sive stars. 



6 DISCUSSION 

The potential for binary systems to merge and form the 
most massive stars is an intriguing solution to the difficulty 
of overcoming the radiation pressure of the forming high- 
mass star. A binary system with a semi-major axis of a few 
AU can be easily perturbed by an encounter with a third star 
that passes near the binary. This encounter radius is signifi- 
cantly larger than the size-scale for direct collisions and thus 
drastically reduces the requirement of an ultradense cluster. 
The encounter time for an interaction at i? C nc is given by 



= 16^/ivnv disp Rl nc ( 1 + o , 



GMbin 



2 W disp^ enc 



(13) 



where w disp is the velocity dispersion and n is the stellar 
density. If we consider an eccentric binary of semi-major 
axis 1 AU, the stellar density n of order 10 6 stars pc~ 3 is 
required to have another star pass within i? cnc « 2 AU, of 
order the apastron separation, in t cnc < 10 6 years, assuming 
fdisp in the range of a few to 10 km s _1 . This is still a 
fairly high density but not drastically so relative to that of 
the cores of young stellar clusters (McCaughrean & Stauffer 
1994), especially considering that they could have been much 
denser in their earlier evolutions (Kroupa, Hurley & Aarseth 
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2001; Scally, McCaughrean & Clarke 2005). Central stellar 
densities an order of magnitude higher are deduced in the 
more massive clusters near the Galactic centre (Figer et al. 
1999). 

An important implication of this work is that the most 
massive stars are more likely to be single stars, or at least not 
in a close binary system. This is due to their being merger 
products and thus the binary system is destroyed in forming 
this very massive star. Intriguingly, well determined masses 
for stars in binary systems rarely exceed w 50 Mq (Bagnuolo 
et al.1992; Rauw et al.1996; Massey et al.2002; Bonanos et 
al.2005) with the most massive stars in a binary being a pair 
of « 70 M (Rauw et al.2004; Bonanos et al.2004). This is 
in contrast with deduced masses for single stars that can ex- 
tend to 150 MQand beyond (Davidson & Humphreys 1997; 
Figer et al.1998; Eikenberry et al.2004). This apparent dis- 
crepancy can be resolved either by the spectroscopic masses 
being far in excess of the dynamical masses, or, by the most 
massive stars having formed from the merger of two massive 
stars in a close binary system. 



7 CONCLUSIONS 

The formation of high-mass close binary systems is a natural 
outcome of forming massive stars through gas accretion in 
a clustered environment. Competitive accretion in clusters 
naturally forms closer, more massive binaries out of lower- 
mass wide systems. These systems can evolve from solar- 
mass binaries with separations of order 100 AU to masses 
of 10 — 3OM with separations of order 1 AU or less. The 
frequency of high-mass binary systems formed in a numerical 
simulation of a forming stellar cluster is 100 per cent for 
stars with m > lOM© and 60 per cent for stars with m > 
5Mq . The high frequency of binary systems and the fact that 
they generally comprise two high-mass stars is in agreement 
with observations (Mason et al.1998; Garcia & Mermilliod 
2001; Preibish et al.1999). The binaries typically have very 
eccentric orbits such that their periastron separations can 
be comparable to, or less than, the stellar radii. Mergers of 
these systems are therefore expected and can be a significant 
factor in forming the most massive stars. Stellar densities 
required to perturb such a binary and force it to merge are 
estimated to be of order 10 6 stars pc~ 3 . 
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